Incompressible Quantum Liquids and New Conservation Laws 
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In this letter we investigate a class of Hamiltonians which, in addition to the usual center-of- 
mass (CM) momentum conservation, also have center-of-mass position conservation. We find that 
regardless of the particle statistics, the energy spectrum is at least q-fold degenerate when the filling 
factor is p/q, where p and q are coprime integers. Interestingly the simplest Hamiltonian respecting 
this type of symmetry encapsulates two prominent examples of novel states of matter, namely the 
fractional quantum Hall liquid and the quantum dimer liquid. We discuss the relevance of this class 
of Hamiltonian to the search for featureless Mott insulators. 



In the twentieth century the basic notions of "symme- 
try" and "order" have set a paradigm that gave rise to 
a lasting cross-fertilization between various branches of 
physics. Landau [1] first put forth the idea that in the 
process of ordering a system also becomes less symmet- 
ric. This phenomenon is commonly known as "symme- 
try breaking", and underlies our general understanding 
of order in a wide variety of systems. Ordering is such a 
ubiquitous tendency of nature that it is difficult to give an 
example where a system remains unordered to the lowest 
temperatures. In recent years condensed matter physi- 
cists have been asking whether the electrons in a solid 
can remain unordered at absolute zero temperature. 

Following Anderson's proposal of the "spin liquid" [2], 
a novel magnetic state with no order, for the past 
twenty years condensed matter physicists have expended 
tremendous effort searching for unordered electronic 
states in solids. [3-5]. Recently there has been signif- 
icant theoretical progress in demonstrating the stabil- 
ity of some examples of such states, and hence prov- 
ing that they are in principle possible [6-8]. In addition, 
it has been shown that there is a strong tie between 
such states and the phenomenon of "quantum number 
fractionalization" . [9] Concurrent with these efforts there 
has been an extensive pursuit for microscopic models ex- 
hibiting these properties. [10-15] The traditional guide- 
line for such pursuits has been "frustration" (see, e.g. 
[16-18]). In this letter we adopt a different viewpoint, 
and investigate whether quantum disordered fermion or 
boson states can arise from unusual conservation laws. 

It is commonly perceived that one of the most difficult 
situations for a system to be free of symmetry break- 
ing is in fractionally filled Mott insulators. For bosons, 
the filling factor is just the average number of particles 
per crystalline unit cell, and for electrons it is half that 
quantity. In this paper, we view any non-disordered in- 
sulating state as a Mott insulator if the insulating prop- 
erty cannot be attributed to the Pauli exclusion princi- 
ple. Whether a fractionally filled Mott insulator can be 
"featureless" (i.e. unordered) has been a central issue of 



debate in recent years. So far all known such systems 
possess order[19]. For example, a half-filled square lat- 
tice of boson atoms is only known to be Mott insulating 
when the atoms localize in a checkerboard pattern and 
hence break the symmetry of the original lattice [20] . 

Physicists have only recently begun to understand 
the reason why symmetry breaking usually accom- 
panies Mott insulators. This progress is due to 
generalizations [21, 22] of an earlier theorem by Lieb, 
Schultz and Mattis.[23] In particular, Oshikawa[21, 24] 
has argued that when a fractionally filled Mott insula- 
tor has an energy gap between its ground state(s) and 
excited states, it must have ground state degeneracy. In 
particular, if the filling factor is p/q, the ground state 
must be at least g-fold degenerate. Empirically, this de- 
generacy is always achieved by breaking symmetry. In 
the boson Mott insulator example given above, a 2-fold 
degeneracy arises because there are two distinct, comple- 
mentary, checkerboard-ordered patterns. For the time 
being we will use the phrase "Mott insulator" to refer to 
a fractionally filled Mott insulator with an energy gap. 

Symmetry breaking is one way to produce the degen- 
eracy required by Oshikawa's result. However if there 
were a way to produce this degeneracy independently of 
symmetry breaking, it would provide a candidate route 
to finding featureless Mott insulators. In the following 
we show that imposing center-of-mass (CM) momentum 
and position conservation is one such way. 

Theorem: Simultaneous conservation of center- 
of-mass momentum and position guarantees q- 
fold degeneracy of the energy spectrum. Consider 
a D-dimensional lattice. We impose periodic boundary 
conditions in a certain direction, say the x-direction, of 
the lattice. Let L be the spatial period in this direction, 
and T be the operator of translation by one lattice con- 
stant in the x-direction. The "cross section" C of the 
lattice is defined such that the total number of unit cells 
is given by CL. Following Oshikawa, let us now consider 
filling factor v — p/q, and C relatively prime to q. In this 
case Oshikawa's argument tells us that the ground state 
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is at least q-fold degenerate if there is an energy gap. 

The (exponentiated) CM position in the x-direction 
modulo L is given by 

r 

The simultaneous conservation of the CM momentum 
and position implies [H, U] = [H, T] = 0. Because of this, 
we can choose the eigenstates of H to be simultaneous 
eigenstates of U. In addition since 

T- 1 UT = cxp 



Here 
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e i 2w C| ^ (2) 



we obtain 
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(3) 

Now we show that q consecutive actions of T on any en- 
ergy eigenstate generates q degenerate orthogonal states 
differentiated by their eigenvalue with respect to U . Let 
the eigenstate in question be denoted by \e lc ^,E > where 
e 1 ^ is the eigenvalue with respect to U and E is the en- 
ergy eigenvalue. Because the translation operator com- 
mutes with the Hamiltonian, the state T\e 1 ^ \E > is also 
an eigenstate of H with the same energy. However its 
eigenvalue with respect to U is 



U T|e 40 , E > = e 427rG i TU\e 1 ^, E > 

= (e^h^ T\e l<t> ,E > . (4) 

The different eigenvalue with respect to U implies that 
the two energy eigenstates \e lc ^,E > and T\e %< ^,E > are 
orthogonal. We can perform this operation q times before 
we arrive at a state with the original U eigenvalue. Hence 
there is at least a 9- fold degeneracy. 

Motivated by the above observations, we would like to 
analyze Hamiltonians featuring CM position and momen- 
tum conservation to determine if they lead to featureless 
Mott insulators. The simplest non-trivial such Hamilto- 
nian is 

H = 9(*,y) C R+X C R _ X C R _ y C R+y . (5) 
R,x,y 

This Hamiltonian hops a pair of particles while con- 
serving their CM position. In general it is quite diffi- 
cult to solve for the eigenstates of this type of Hamilto- 
nian. However, for a special one-dimensional case some 
insight can be gained due to a connection with a well- 
studied problem in physics, the fractional quantum Hall 
effect. To be specific, we consider the the Trugman- 
Kivelson (TK) Hamiltonian[25], which has Laughlin's 
1/3 wavefunction[26] as the ground state. When studied 
on a torus with dimensions L x and L y , the TK Hamilto- 
nian can be written in terms of the creation/annihilation 
operators of the lowest Landau level as follows [27]: 



H 



E 

R.x.y 



f*{x)f{y)C R+x C R _ x C R . y C R+y , 
f(x) = k 3/2 Y,( x - nL)e- K ^ x ' nL ^ 2 . (6) 



L = L x Ly/2nl 2 3 , k = 2ttIb/ L y . 



(7) 



In the above lg = yJhc/eB is the magnetic length. This 
Hamiltonian acts upon a system of spinless fermions on 
a ring of L sites at filling factor 1/3. 

One important consequence of the connection between 
Eq. (6) and the quantum Hall effect on a torus is an 
unusual property of Eq. (6), which we call duality. By 
performing Fourier transforms it is simple to show that 
when written in terms of operators that create and anni- 
hilate particles in a fixed momentum state, Eq. (6) reads 

H = Y, f*( k )f(l) Cg+fe^Q-fe C Q - q C Q+q , (8) 

Q,k,q 

where f(k) = J2 X e i2kx f{x). For the / used in Eq. (6) / 
has the same form as /. This implies that the Hamilto- 
nian in the momentum space is also described by Eq. (6) 
except k — > jjj. The factor 2ir/L arises from the lattice 
constant in the reciprocal space. For each L there is a 
special k value (k* = y/2n/L) for which the real space 
and momentum space k are the same. The duality im- 
plies that energy spectra at k and 2tt/kL are identical. 
From the perspective of the quantum Hall liquid on a 
torus, the duality merely signifies that interchanging L x 
and L y leads to the same physical system. 

Analyzing the Hamiltonian Eq. (6) we find that for 
all practical purposes the ground state is a featureless 
Mott insulator when k is small. In addition, we find the 
surprising fact that on a torus with finite circumference 
in one direction (but infinite circumference in the other) 
the Laughlin liquid has non-vanishing density wave order, 
and is adiabatically connected to a Wigner crystal. 

First we focus on the large k limit. In the large k limit 
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we can expand the Hamiltonian in the parameter e~ K . 
As a result Eq. (6) is accurately approximated by 



H 



L 

E 

i=l 



[/(l/2) ; 



rii+xrii + f(l) 2 n i+2 n t 



(9) 



where ni = C^Ci. This Hamiltonian imposes energy 
penalties for having nearest neighbor and next-nearest 
neighbor particles. Its ground states are "Wigner crys- 
tals" where one out of every three lattice sites is occu- 
pied. Clearly at v = 1/3 there are three such ground 
states. These ground states possess density wave order 
(hence are not featureless). 

To quantify the degree of crystalline order we introduce 
the order parameter O associated with a spatial period 
of three lattice constants: 



1 L 

o = -Y 



< nj >. 



(10) 
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In the above N is the total particle number, and < rij > is 
the expectation value of the site occupation in (any one 
of) the ground state(s). O is normalized such that its 
modulus becomes unity in the extreme crystalline limit. 
In Fig. (l)a) we plot O as a function of re, evaluated nu- 
merically for system sizes L — 6, 9, 24. For re > 1.5 we 
obtain O ~ I implying that the ground state is nearly 
a perfect crystal. As expected, as re decreases the or- 
der parameter decreases and practically vanishes below 
re rs 0.5. [28] This finding raises an important question: 
Is the regime below n = .5a featureless Mott insulator? 
The weak size dependence of the crystalline order pa- 
rameter near re = 0.5 is inconsistent with the existence 
of a continuous phase transition, which is required if the 
state below re = 0.5 is truly free of crystalline order. The 
smooth evolution of O with re rules out the possibility 
of a first-order phase transition. Indeed, while the order 
parameter becomes exponentially small at small re (nu- 
merical precision limits our study to re > .25.), a careful 
analysis of the numerical data indicates that at any finite 
re, it will not vanish as L — > oo. Using standard methods 
of extrapolation [29] we find that the order parameter for 
L — > oo is given by exp(— l/r(re) 2 ), where r(re) vanishes 
linearly at re = (see inset in Fig. (l)(a)). According to 
this result, density wave order exists as long as re ^ 0. 

In terms of the quantum Hall connection the fact that 
the Laughlin liquid on a torus with any finite L y (and 
infinite L x ) has a nonzero density wave order is very sur- 
prising (to us at least). From the perspective of a ID 
lattice problem, however, it can be shown that in general 
there exists a local order parameter that can be used to 
distinguish the degenerate ground states required by Os- 
hikawa's argument. [30, 31] However, e.g., for re — 0.25 
the crystal order is 15 orders of magnitude weaker than 
that at re = 1.5! Hence for all practical purposes this is 
a state without order. 

In the following we carefully study the energy gap as 
a function of re to demonstrate that 1) the system is a 
Mott insulator for all re, and 2) there is no quantum phase 
transition between the large and small re regimes. In 
Fig. (l)b) we show the numerically computed energy gap 
vs. re for system sizes up to L = 30. Only data for 
even particle number (to avoid even/odd effect) are used. 
The limit L — ► oo and re fixed amounts to studying a 
Hall torus with L x — > oo and fixed L y fixed. To avoid 
finite size effect from small L x we restrict L x > 121b 
or equivalently, kL > 12 (for L — 30, re > 0.4). With 
these restrictions the data points collapse onto a single 
curve, as shown in Fig. (l)(b). Note that according to 
this result, the gap remains robust even in the "transition 
regime" around re = 0.5. According to quantum Hall 
physics the energy gap should not depend on L x or L y 
when both are much greater than lg. Thus for re < 0.4 (or 
L y > 5irls for the Hall torus), we expect little variation 
in the extrapolated gap value. The asterisk in Fig. (l)(b) 
is obtained by extrapolating the gap value at k*(L) to 




FIG. 1: The density wave order parameter and the energy 
gap of Eq. (6) as a function of re. Data are obtained for 
L < 30. a) The crystalline order parameter. Inset: The 
quantity Y = (— log 0) -1//2 , extrapolated to L = oo using 
the alternating epsilon algorithm. [29] The smooth curve is a 
polynomial fit of the data for .25 < k < .45 (crosses), b) The 
energy gap for data restricted to kL > 12 and even particle 
number (see text). Inset: The energy gap is plotted versus 
Of{\) 2 . Blue and red symbols are for re > 1.1 and re < 1.1 
respectively. Two distinct branches are clearly visible. 



L = oo (note that liniL^oo K*(L) = 0). 

From these results the energy gap is nonzero for all 
k.[32] However, in the inset of Fig. (l)b) we demonstrate 
that two different mechanisms cause the energy gap at 
small and large re. For large re the gap is the energy 
penalty of having next-nearest-neighbor particles. The 
latter quantity is proportional to f(l) 2 0. However the 
energy gap plotted against f(l) 2 exhibits two branches, 
indicating a different mechanism for the energy gap at 
small re values. Hence the small re state is not just an or- 
dinary crystal with extremely small order. Thus even in 
the most unfavorable case (ID) , our approach has suc- 
ceeded in producing a Mott insulator that is practically 
featureless. In higher dimensions, where there are no 
theorems forbidding featureless states, we expect Eq. (5) 
will lead to truly featureless Mott insulators (see later 
discussion on quantum dimer models). 

Our analysis also leads to important implications for 
the fractional quantum Hall liquid on a torus. In the 
quantum Hall community, it is well known that when L y 
is large compared to lg (i.e. when re is small), the Laugh- 
lin state on the torus should be practically indistinguish- 
able from that in the infinite 2D plane (L — > oo, re = 0). 
Quite unexpectedly, however, our findings imply that the 
Laughlin state in this regime is adiabatically connected to 
a state with strong crystalline order at small L y .[33] The 
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above result reveals a danger in over-interpreting adi- 
abatic continuity, because the Laughlin state is clearly 
qualitatively different from the electron crystal. How- 
ever, since discrete quantum numbers are preserved by 
adiabatic continuity, both the ground state degeneracy 
and the quasiparticle charge are the same in the quan- 
tum Hall and the crystal state. The obvious threefold de- 
generacy of the electron crystal originates from the three 
distinct center-of-mass positions of the system. The ex- 
cited states of the crystal are domain walls carrying ±1/3 
charge due to the Su-Schrieffer counting argument [34] . 

Our findings are consistent with the notion that the 
Laughlin state on a finite torus has topological order. As 
proposed by Niu and Wen[3], unlike in a symmetry bro- 
ken state, disorder lifts the ground state degeneracy of 
a topologically ordered state by an amount that is ex- 
ponentially small in the system size. In our case, such 
exponential dependence comes from the weak order pa- 
rameter itself. Based on this we conclude that despite the 
common perception, the presence of symmetry breaking 
order should not be taken as excluding topological order. 

Our belief that Eq. (5) can lead to truly featureless 
Mott insulators in higher dimensions is supported by the 
following connection to the quantum dimcr model. [10] 
The quantum dimcr model is a special case of the bosonic 
version of Eq. (5) in two dimensions. Indeed, once dimcrs 
are reinterpreted as point bosons residing at bond cen- 
ters, the quantum dimcr Hamiltonian becomes a center- 
of-mass position conserving model described by Eq. (5). 
Moessner and Sondhi[ll] argued that the quantum dimcr 
model on a triangular lattice exhibits an unordered quan- 
tum phase. This supports our belief that Eq. (5) leads 
to truly featureless Mott states in higher dimensions. 

Finally, although we have focused on Mott insulators, 
gaplcss states described by Eq. (5) are equally interest- 
ing. We can prove that such a gapless system cannot be 
an ordinary metal (fcrmions) or a supcrfluid (bosons), 
because center-of-mass position conservation implies the 
absence of Drude weight /superfluid density. For example 
in one space dimension such a gapless liquid will not fall 
within the usual Luttinger liquid paradigm. 

To conclude we have presented the idea that simulta- 
neous conservation of CM position and momentum can 
lead to fractionally filled featureless Mott insulators. It is 
remarkable that the simple model given by Eq. (5) uni- 
fies systems as diverse as the quantum Hall liquid and 
the quantum dimer liquid, whose effective field theories 
are as different as Chern-Simons gauge theory and Z2 
gauge theory, respectively. We hope our results will spur 
the exploration of new directions and lead to a wealth 
of new states of matter that have, so far, escaped our 
attention. 
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